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Abstract 

We use two non-perturbative methods to obtain the anisotropy derivatives of 
the coupling constants (the anisotropy coefficients) of 577(3) lattice gauge theory. 
These coefficients appear in the derivative formulae for the energy density and the 
pressure. We calculate them for the standard Wilson and two improved actions, 
the (2 x 2) and the Square Symanzik action. Both methods lead for all investigated 
actions to compatible results which are clearly different from their known asymptotic 
perturbative limits. With increasing f3 the limits are however approached in all cases. 
Our findings strongly support the equivalence of the integral and derivative methods 
for the calculation of energy density and pressure. 



*Present address: IBM, P.O. Box 265, 00101 Helsinki, Finland. 



1 Introduction 



For the investigation of the phase transition from ordinary hadronic matter to the 
quark gluon plasma and the approach of the plasma to its high temperature limit 
the calculation of the energy density e and the pressure p is of central importance. 
These observables are defined as derivatives of the partition function Z with respect 
to temperature T and volume V 

1 dlnZ dlnZ 

P = T ^77- ■ (1-1) 



Vd(l/T) ' r dV 

In the lattice formulation of QCD the temperature and the volume therefore have 
to be independent variables. On a lattice with a fixed number of points N a in 
each spatial direction, and N T points in the temporal direction, this is achieved by 
choosing corresponding independent lattice spacings a a and a T so that 

V = (N a a a ) 3 and T' 1 = N T a T . (1.2) 

As a consequence two different coupling constants have to be introduced for the 
spatial and temporal contributions to the action. For SU(N C ) lattice gauge theory 
we consider actions of the following form 

s = p„ E s^W+Pr E > (i-3) 

X,I1<1><4: X,fJ,<4 

where S^ u {x) is a generalized plaquette in the (/i, v) —plane. Generalized plaquette 
expectation values are then defined by 

c t x,n<u<A 

P r = E *W*)> • (1-5) 
oiv <T iv -r x,fj,<4 

The couplings /3 CT)T may be rewritten in terms of the coupling constants g 2 ar as 

Pa = —5- • 7 and (3 T = — =- • f . (1.6) 

Here ^ = a a /a T is the anisotropy. One may also express the couplings with the bare 
anisotropy 7 and a common coupling /3 through 111] 



or equivalently 



Pa = — , Pt= Pi 

1 



X 



The couplings depend on a a and a T , or, alternatively, on a a and £. Following ref. @ 
we adopt the latter two as independent variables. The energy density and pressure 
are then 

-fM' f - 3p= v( a '(K^ (L9) 

Obviously, the energy density contains now terms, which are proportional to deriva- 
tives with respect to the anisotropy £, because 



dPr 



-Pa + 2N ( 



d 9 g 2 

p T +2N c e dg;2 



1.10) 

1.11) 



<9£ <9£ • 

After taking the derivatives one can evaluate the expressions on isotropic lattices 
At £ = 1 we have 

7 = 1 > 9a = 9r = 9 , a a = a T = a , (1.12] 
and the relation for the (3— function [0 



da 



;i.i3) 



Energy density and pressure are now conveniently combined in the two equations 



e + p 
e — 3p 



8N c N*g 



4„-2 



1 _ f^k! _ ^ 



e=i. 



{Pff - Pr) 



12N C N* ( %! + [2P - (P CT + P T )] . 



<9£ <9£ 



;i.i4) 
;i.i5) 



Here Pq denotes a generalized plaquette on a symmetric (T = 0) lattice. Its inclusion 
normalizes both e and p to zero at T = 0. 



We define the anisotropy coefficients follows 



;i.l6) 



In the weak coupling limit, they are related to the derivatives of the anisotropic 
parts of the couplings 



g£(a,t)=g-\a) + c»Jt)+0(si 1 



;i.i7) 



by 

The quantities (£) were calculated in ref. for the Wilson action perturbatively. 
For the difference one obtains 

N 2 - 1 

c ff (0) - c r (0) = 0.146711 - N c ■ 0.019228 , (1.19) 

l \ r 



and, from the renormalization group equation and eq. (|1.13 ) 
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c CT (0) + c T (0) = 6 = i ^|. (1-20) 

These perturbatively calculated anisotropy coefficients were used in early finite tem- 
perature lattice calculations. In SU (3) gauge theory, however, this led to unphysical 
behaviour of the pressure: it could become negative and at the critical point a gap 
in p was observed |4], [5| . A second, rather important consequence of taking only the 
perturbative coefficients is the rapid approach of e and p to their respective ideal gas 
limits soon after the transition, though the limits themselves are independent of c CT)T . 
The actual form of the approach is of course dependent on the residual interaction 
in the quark gluon plasma, and therefore of physical relevance. 

In ref. an alternative method for the calculation of energy density and pressure 
was proposed, which is now in common use. Here the pressure is obtained without 
the anisotropy coefficients from an integral over plaquettes. In order to be consistent, 
the integral method should be equivalent to the use of non-perturbative anisotropy 
coefficients, i.e. the corresponding values for c CT)T should be essentially independent 
on finite size and cut-off effects. First results in ref.Q showed however considerable 
cut-off effects for the Wilson action. To clarify the situation, we started || a non- 
perturbative calculation of the coefficients with the matching method, which we 
want to complete in this paper. 

In the following section we discuss the integral method in some more detail as 
in and || and apply it to the Wilson action. The matching method is described 
and applied to the same action in section 3. For comparison we determine in section 
4 the anisotropy coefficients as well for two improved actions, the (2 x 2) and the 
Square Symanzik action M. We close with a summary and the conclusions. 
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2 The integral method 



For large homogeneous systems with isotropic interactions, the volume derivative in 
the pressure formula (|1.1|) may be replaced by a volume division. The pressure and 
the free energy density are then related by 

p=-f = Ll nZ . (2.1) 
Instead of calculating InZ directly, one may take advantage of the relation 

d{3 W' 1 J 

and integrate the measured expectation value over a (3— range. For £ = 1 we obtain 



then 

P 

p p 

rp4 ~ rp^ 



+ J dp'[2P - (P CT + P T )] . (2.3) 



* ft 

As in eqs. ( |1.14|) and Ql.151) , the inclusion of the P — term serves to normalize p at 
T = to zero. Obviously, the last equation will lead to a continuous pressure, even 
if there is a gap in the plaquettes at the critical point. The energy density is found 
from 



e - 3p ^ Ar . dg 



-2 



, r -Q Nc N*a^[2P - (P a + P T )\ , (2.4) 

by just adding 3p/T 4 . Here we have replaced the sum of the coefficients by the 
(3— function. This function and/or a = a(g 2 ) have to be known non-perturbatively, 
not only to evaluate the last equation, but also to determine the temperature T = 
l/N T a on fixed N T — lattices at any given coupling. Parametrizations of the (3— func- 
tion for the SU(3) Wilson action have been derived by several groups JF], [U], |TT 



They are rather similar above (3 = 5.9, but below this value they differ from each 
other. Probably this is due to the transition from the weak coupling to the strong 
coupling regime. 



2.1 The high temperature limit 

The two methods for the calculation of e and p on the lattice suffer from cut-off 
effects, especially for small N T — values. Though these cut-off effects are of similar 
size, there is a non-negligible difference. For the same method, the high temperature 



4 



N T 


R D 


Ri 


r = R D /Rj 


4 


1.4952 


1.3778 


1.0852 


6 


1.1816 


1.1323 


1.0435 


8 


1.0867 


1.0659 


1.0195 



Table 1: The ratios R(N T ) = e(N T )/esB = p(N t )/psb for the Wilson action calcu- 
lated with the derivative (Rd) and the integral methods (Rj). 

limits of e and p have the same cut-off dependencies, because in both cases the same 
formula for (e — 3p)/T 4 is used and 

(e - 3p)/T 4 = 0(g 4 ) . (2.5) 

However, as different methods are based on the evaluation of different operators, 
which have different cut-off dependencies, the high— T limit derived in different 
approaches leads also to distinct cut-off dependencies. In the derivative method the 
high temperature limit is derived from the lowest order of the g 2 — expansion of the 
energy density (£ = 1) 

^ = W c N*g- 2 (P a -P T ) + 0(g 2 ) , (2.6) 
= l(N 2 - l)iV r 4 (p( 2 ) - PW) + 0(g 2 ) . (2.7) 

The quantities Pjf} are the corresponding expansion coefficients of the generalized 
plaquette expectation values 

P«,r = g 2 ^^P%+0(g 2 ). (2.8) 
For the Wilson action one obtains the following deviation of the energy density, 



e(iV T ), in the high T limit from the continuum result, e^^ HI 



If the integral method is used, the calculation of the high temperature limit starts 
from eq. (|2.1|) and the lowest order term of InZ. Here the T = subtraction has to 
be done explicitly. With the Wilson action one arrives at 



where of course €sb = —3fsB- The full 1/N 2 — expansions can be calculated nu- 
merically. In Table p] we list them for the two methods in the case of the Wilson 
action and N T = 4, 6, 8 . We show also the ratio r between the two high temperature 
limits. We stress that this difference, which expresses itself in the ratio r, will play a 
crucial role in the discussion of the anisotropy coefficients presented in the following 
sections. 



2.2 Determination of the anisotropy coefficients 



As mentioned already in the introduction, the integral method implies the use of 
non-perturbative anisotropy coefficients. We can determine them by comparing the 
pressure formulae from the two different approaches 



T*J D 



p_ 

T 4 / / 



(2.11] 



and by exploiting the (3— function relation, eq.( |1.13|) . In the pressure equation (|2.11 



we have inserted the factor r as a compensation for the differing cut-off effects. In 
principle, r = r(g 2 ) . We approximate r by its value at g 2 = (Table |I]). From 
eqs.( |1.14"|) and ( |1.15| ) we obtain 



p_ 

T 4 / z > 



N r Nt 



(P CT - P T ) - 3( Cff + Or) (2P - (P CT + P T 



[2.12) 



and therefore 



a dg 



4 da 



-c a - 



Pr 

a dg 



P 



2NM VT 4 



3a do 2 / 

f + Tir( 2P »-( p " + ^ 



2 da 



(2.13) 
(2.14) 



In the following we discuss the behaviour of c a . The behaviour of c T is similar, 
because both coeffcients are related by the smooth j3— function. It is clear, that the 
numerical determination of the coefficients becomes problematic below the transi- 
tion point, because there both the pressure and the plaquette differences are small. 
Indeed, this is seen in Fig. f[[ where we show c a for the SU(3) Wilson action as ob- 
tained from eq. fl2.13| ) on N T = 4, 6, 8— lattices using interpolations of the \P a — P T ] 
and [Pq — (Pa + P T )/2]— data of ref. 0. Apart from the close vicinity of the critical 
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Figure 1: The coefficient c CT for the SU(3) Wilson action from plaquette interpola- 
tions. The dashed curve is without cut-off corrections from ref. ||. Also shown is 
the perturbative result for (3 — > oo , c a ^ and the sum c a + c T . 
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Figure 2: The coefficient c a for the S'L r (3) Wilson action with errors coming from the 
plaquette data . The solid curve is the fit, eq. ( |2.i5|) , the yellow line a corresponding 
fit to data (not shown), which one obtains if p/T 4 (/3 = 5.65, N T = 4) = 0.05 . 
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(3— values the results for different iV T are consistent with each other. For comparison 
we show also the previous result for N T = 4 from ref. f|, where the difference in 
the cut-off effects of the two methods was not taken into account, that is r = 1 was 
used. Obviously the influence of the r— factor is important. Moreover it leads to 
slightly lower results for c CT , though, as before, there remains a substantial deviation 
from the asymptotic value c a>00 = c CT (0) from eqs. ( |f .f 9j ) and ( |f ,20| ). 



It is difficult to estimate all possible error sources in the calculation of c CT , as 
for example that of the parametrization of the (3— function. The major errors are 
certainly due to the errors in the plaquette expectation values. In Fig. |2] we present 
the results including these errors, but omitting the respective critical regions. We 
can summarize these data with a simple Pade fit of the form 

I + diq 2 + d 2 q i 

— , (2.15) 



where d = -0.64907, d 1 = -0.61630, d 2 = 0.16965 . 

A further source of errors is the result for (p/T 4 )/ . Whereas the integration itself 
contributes only a negligible error, the size of the pressure at the lower integration 
boundary (3q is unknown. We have identified f3$ with the point where the integrand 
[Po — (P CT + -P T )/2] is disappearing and assumed, as usual, that there p/T 4 is zero. 
This may be not true. In order to assess the effect of a non-zero pressure contribution 
we have repeated our calculation with a test value p/T A = 0.05 at (3$ = 5.65, T /T c = 
0.917 for N T = 4. As a result the value of c CT is diminished by 0.022 in the whole 
(3— range considered. The same shift in c a is observed also for N T = 6 and 8. For 
comparison we have plotted in Fig. ^| a fit of the form ( |2.15| ) to the shifted data as 



well. Here we find d = -0.60983, di = -0.78889, d 2 = 0.36164 . 

Some insight into the question of the pressure size below the transition is gained 
from the behaviours of the plaquette differences around T c in 577(3) and 577(2). In 
Fig. H we compare them to each other. The main difference between the two gauge 
groups is evidently arising from the fact, that the respective phase transitions are of 
different order. Also, there is no doubt, that in 577(2) the pressure is non-zero well 
below T c . 

As we have seen in 577(3) the inclusion of the cut-off factor r considerably lowers 



the N T = 4— result for c CT . In SU(2) the values which were previously [14] determined 
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Figure 3: The plaquette differences in the vicinity of the deconfinement transition 
for SU(2) and ST/ (3) on N T = 4-lattices. The S77(2)-data are from ref. O. 
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Figure 4: The coefficient c a for the SU(2) Wilson action with errors coming from 
the plaquette data and r = 1.0852. The dashed curve (r = 1) and the sum c a + c T 
are from ref. [14[], the red points from ref. Also shown is c^oo . 



for c CT with the integral method and r = 1 are well in accord with points measured 



in an independent way at (3 = 2.30 and (3 = 2.375 by Ejiri et al. |15 |. One has to 
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check then, whether this agreement is lost if the correct r is taken into account. We 
do this by using the new high precision SU(2) data of ref. [|L3|1 , already shown in 
Fig. ^ the (3— function of ref. [14] and the r— value from Table |l|. The resulting 
Co- is plotted in Fig. f|. We find only a small change in this (3— range, the result 
is lowered at (3 = 2.30 by 0.0067 and at (3 = 2.35 by 0.0142 . The reason for this 
is, that up to (3 — 2.35 the difference [Pq — (P a + P r )/2] is still large and only at 
higher couplings the pressure term with the factor r is dominating in eq. (|2.13|) . 
We expect therefore a stronger decrease at higher (3— values. 



3 The matching method 

This method was first used by Burgers et al. [ff] to measure the anisotropy £ as a 
function of the bare anisotropy 7 on anisotropic lattices. With this information one 
may as well determine the anisotropy coefficients. Consider the couplings 

9a 2iV c 7 ' 9r 2N c £ ' 1 ' ' 

and their derivatives with respect to the anisotropy £ 

d 9 - 2 _ _ a / 1 Sin/3 01n 7 \ 

d9r 2 2 I 1 51n/3 «91n 7 \ ,„ 



At £ = 1 we may apply eq. (|1.12j ) to find the difference 

$7 



c a -c T = 2g- 2 [1 - ] , (3.4) 



and, using eq. (|1.13|) , we obtain 



_,, 1 i 6^7^ adg z 



^ = ±e ^-sej^-sir- (35) 

A measurement of the function £(7) in the neighbourhood of £ = 1 will therefore 
enable us to calculate the anisotropy coefficients, once we know the (3— function. In 
principle, the anisotropy £ may be determined with two measurements of a physical 
observable which depends on a distance, which can be chosen in a spatial or in the 
temporal direction. At the same physical distance the expectation values have to 
be the same and thereby fix the anisotropy. This is the idea of matching. 
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3.1 Matching of Wilson loop ratios 



Suitable quantities for the matching process are obtained from Wilson loops of size 
nx x n 2 (in lattice units) . The Wilson loops are related to the heavy quark potential 

W(xi,x 2 ) ~ exp[-xxVi(x 2 )] , (3.6) 

for x\ — > oo . Here, Xi = n^ai. The potential V\ differs from the continuum potential 
V by a term Vq, the self-energy of the heavy quarks, which is dependent on lattice 
spacing, but independent of x 2 0]. The natural way to proceed then is to build 
ratios of Wilson loops, which depend only on Vi(x 2 ) for large x\. We use the following 
ratios 

Ri{ni,n 2 ) = — — — ; 3.7 

W(x 1 ,x 2 ) 

W(xi + 2a 1 ,x 2 )W(x 1 + ai,x 2 ) 

R ^ n ^ = • (3 - 8) 

On an anisotropic lattice there are two different types of loops : space-space {W aa ) 
and space-time (W aT ) Wilson loops. The corresponding ratios R a and R T are mea- 
suring the same potential V(x 2 ), if the matching condition for the physical distance 

x 2 = y = t or n t = £n y , (3.9) 

is fulfilled, apart from an n— independent factor k = k aa /k aT , which is due to the 
dependence of Vb on a a and £ 

Ra(n x ,n y ) = k ■ R T (n x ,n t = £n y ) . (3.10) 

At £ = 1 the factor k is of course 1 . We do not implement any smearing for 
our Wilson loops, because then space and time links would have to be smeared by 
the same amount in physical units. However, in order to enhance the accuray of 
the Wilson loop measurements we are using the fast link integration technique of de 
Forcrand and Roiesnel [0 wherever this is feasible. It cannot be done for loops with 



an Hi = 1 and at loop edges. In applying this algorithm we are saving in addition 
half of the computer time by taking advantage of the fact that the necessary link 
contour integrals are real and can as well be obtained from a half contour in the 
complex plane. 
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For the Wilson action simulations were performed at four j3— values: 5.4, 5.7, 6.3 
and 7.2 using a 16 4 lattice for 7 e [0.92,1.08], a 16 3 x 32 lattice for 7 e [1.1,2.0] 
and a 16 3 x 48 lattice for 7 = 3.0. The Wilson loops were calculated after every 
fourth sweep through the lattice, where a sweep consisted of one heatbath update 
succeeded by four overrelaxation steps. The integrated autocorrelation times which 
one finds then are increasing with increasing f3 from 0.5, 1.5, 3.0 to over 3.0 for 
(3 = 7.2 . On the average we took 5000-7000 measurements for j3 = 5.4 — 6.3 and 
600 for f3 = 7.2 . Generally, for the same number of measurements, the relative 
error of a Wilson loop increases with the size of its area. The increase itself is 
much steeper at lower (3— values than at higher ones. Larger area Wilson loops are 
therefore only included in the matching procedure for the higher j3— values. In Fig. 
|5] we show a typical example of this behaviour. The difference in accuracy between 
link integrated and plain loops is clearly seen in the plot. 
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Figure 5: The relative errors of spatial Wilson loops as a function of their 

7 = 1.06. Open symbols show the errors of link integrated loops, filled symbols 

those of not integrated ones. 



After measuring the Wilson loops at a fixed value of (3 and 7, we compute 
the ratios R a and R T . For 7 > 1 we connect the timelike ratios R T with spline 
interpolations. In order to improve the interpolations we include as well the ratios 
with n t = 1 . The spatial ratios R a are then shifted in n y by a factor £ and in height 
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Figure 6: The ratios Ri tT (n x , n t ) at (3 — 6.3 and 7 = 1.5 for fixed n x = 2, 7 plotted 
vs. 2n x + n t and connected by splines. The ratios Ri )a {n X) n y ) are once plotted vs. 
n y and also shifted, that is vs. 2n x + £n y , with £ = 1.63 . 
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Figure 7: The ratios R2, T (n x , n t ) at /3 = 6.3 and 7 = 1.5 for fixed n x = 3, 6 plotted 
vs. 2n x + n t and connected by splines. The ratios R2,a(n x ,n y ) are once plotted vs. 
n y and also shifted, that is vs. 2n x + £n y , with £ = 1.63 . 
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Figure 8: The result for £ from the matching of R\ at /3 — 5.7 and 7 = 1.5 as a 
function of A = m x m y . In the matching all ratios with n x > m x and n y > m y were 
taken into account. The lines show the final result with an error band. 



by a factor k such that the sum of the squared deviation (x 2 ) of the shifted R a — values 
from the respective R T — interpolations becomes minimal. For 7 < 1 the roles of R a 
and R T are interchanged. The fitting is done for all ratios with n x > m x and n y > m y 
at the same time, with suitably chosen minimal values m x and m y . In so far our 



procedure differs from the matching prescription of Klassen [IS], who matches single 
ratios and looks for a possible plateau of shift values. In Figs. |6| and we show 
matching examples at (3 = 6.3 and 7 = 1.5 for the ratios R\ and R2 ■ The matching 
in Fig. H is optimized for m x = m y = 2 , in Fig. |7] for m x = m y = 3 . Both 
lead to the same value £ = 1.63 for the anisotropy. For Ri we find k = 1.01 and 
k = 1.02 for i?2 , that is k is always very close to one. One may as well perform fits, 
as Klassen [0 does, with fixed k = 1. The corresponding \ 2 is then considerably 
larger, the value for £ increases slightly. We shall come back to this point again. 

As has been mentioned already, we have included only ratios in the matching 
process which are built from Wilson loops with at least minimal extensions m x and 
m y . This has been done for two reasons. Once, x\ should anyhow not be too small 
(see eq. ( |3.6| )) and second, Wilson loops with a single link on one side cannot be 
link integrated and are therefore less accurate. That is why we disregard ratios with 
n x = 1 and/or n y = 1 in the final analysis. The influence of the chosen minimal 
values on the matching result for £ is demonstrated in Fig. [| for R\ at (5 — 5.7 and 
7 = 1.5 . We observe, that apart from the m x = 1 data all other measurements for 



14 



£ are consistent with each other. 



3.2 Results for c a<T from matching 

In order to obtain the anisotropy coefficients we measure now the function £(7) in 
the neighbourhood of 7 = 1 with the aim of determining the derivative {d^/dj)^=i . 
In principle one would choose 7 rather close to one. However, the difference between 
spatial and temporal Wilson loops becomes then very small as well and the matching 
will be inexact. Fortunately, it turns out, that £ is linear in 7 in a wide range around 
7 = 1 at all (3— values. A successful strategy is therefore to measure the Wilson loops 
with high precision at not too many 7— values to determine the slope of £(7) . In 
Fig. |^ we show as an example the function £(7) for the Wilson action at (5 — 5.7 as 
obtained from the matching of Ri with m x = 2,m y = 3, together with a linear fit 
to the data points. 



3 . 



2 . 



1 . 





1 1 1 1 


ill. 




1 | 1 1 1 


1 1 






17 1 . 1 

1.0 


| 1 1 1 


| 1 1 1 1 








0.9 


1 1 1 1 
















. 9 


1 . 




1 . 1 

1 



1.0 1.5 2.0 2.5 3.07 



Figure 9: The anisotropy £(7) at (3 = 5.7 as determined from the ratio R\ with 
m x = 2, m y = 3. The red lines are from a linear fit in 7, the blue lines show the 
error band . 

As in the last section, we use now the (3— function of ref. [0 to deduce the 
anisotropy coefficients from the measured derivative (d^/dj)^=x . The corresponding 
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Figure 10: The anisotropy coefficient c a of the Wilson action from the matching of 
R\ (black) and R2 (blue points) as a function of A = m x m y for (3 = 5.4, 5.7, 6.3, 
and 7.2 . Points with m x = 2, 3, 4, 5 are denoted by A, x, y, + , respectively. 
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p 


5.4 


5.7 


6.3 


7.2 






0.413(23) 


0.315(09) 


0.261(08) 


0.227(15) 


C T 




-0.374(23) 


-0.276(09) 


-0.202(08) 


-0.159(15) 


c a {k = 


1) 


0.397(17) 


0.348(06) 


0.296(06) 


0.261(12) 


c T {k = 


1) 


-0.358(17) 


-0.309(06) 


-0.237(06) 


-0.193(12) 



Table 2: The matching results for c a and c r for the Wilson action . 



results for c a are shown in Fig. [10] as a function of the minimal area A = m x m y of 



the Wilson loops, which were included in the ratio-matching. The data for (3 = 5.4 
are somewhat problematic in several respects and can only serve as an indication 
for the true value. As demonstrated already in Fig. [5] too large Wilson loops cannot 
be included in the matching for (3 = 5.4, because their errors are too large. The 
second handicap is the unknown or ambiguous (3— function in the strong coupling 
region. We have therefore just taken the same value for c a + c T at (3 = 5.4 and 5.7 . 
At (3 > 5.7 we find that the results are essentially independent on the minimal area, 
at /3 — 7.2 the errors are larger because of the lower statistics of the measurements. 
The final results for c a and c T are given in Table |2| . We obtained them by fitting the 
single values to a constant with the x 2 — method. Their average errors were estimated 
from the error estimate of the constant times the square root of the number N of 
fitpoints. The error estimate was rescaled whenever x 2 /{N — 1) > 2 . This was only 
necessary for f3 = 5.4 . In order to be able to compare our results to those of Klassen 
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Figure 11: Comparison of the anisotropy coefficient c a of the Wilson action from 
the matching of R\ with k ^ 1 (black) and k — 1 (red points) as a function of 
A = m x m y for (3 — 6.3 . Points with m x = 2,3,4,5 are denoted by A, x, y, + , 
respectively. 
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T8| , we have repeated the matching analysis under the assumption k = 1 . Though 
the single matching results seem to have smaller errors then, a constant x 2 — fit as 
before yields always a \ 2 P er degree of freedom which is larger than 2. The inclusion 
of small Wilson loops in the matching seems to be more crucial here and also to 
lead to higher values for c a . This is observed in Fig. |11], where we compare the two 
methods at (3 = 6.3 . 
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Figure 12: Comparison of all available data for the anisotropy coefficient c a of the 
Wilson action. Further details are explained in the text. 



Finally, we have gathered in Fig. ^ all available results for the Wilson action 
c CT : those which we obtained from the integral method are represented by the fit, 
eq.(|2.15|), and the fit showing the possible influence of p/T 4 (/?o) ; our matching re- 



sults for k 1 and k = 1; the 577(3) data of Ejiri et al. ||15|| and those of Klassen p8 



both for our f3— function. We observe consistent behaviour of the integral and match- 
ing results for k = 1 , those for k ^ 1 are still within the error bars (see Fig. 0). 
There is also full agreement with the data of Ejiri et al. . Only the results of Klassen 
are definitely higher and insofar incompatible with all other measurements of c a . 
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4 Anisotropy coefficients for improved actions 



Up to now we have discussed in detail the calculation of the anisotropy coefficients 
for the Wilson action only. In the following we consider the more general family of 
actions defined by Wilson loops of size up to 4 in the (fi, i/)-plane 



•V = °i,i ( x _ ^T Re Tr :. 

+ ai , 2 (l-^Re^ (a + 



«, 2 (l-iReTrr j.il.l; 



Here, the parameters an, ai 2 and 022 are constrained by the equation 

an + 4a 12 + 16a 22 = 1 , (4.2) 

which ensures the correct continuum limit. The action is 0(a 2 ) improved, if addi- 
tionally 

a 12 = ---, a 22 = -^ u --J • (4.3) 

Still, one parameter, for example an, is free. For an = 5/3 we obtain the Symanzik 
action (a 12 = —1/16, a 22 = 0), for an = 4/3 the (2 x 2)— action (a i2 = 0, a 22 = 
— 1/48). Alternatively, one may require that the propagator can be diagonalized. 
This facilitates the analytic calculation of the anisotropy coefficients in the weak 
coupling limit. In that case we have the condition 

an = / , o 12 = 2za n , a 22 = z 2 a u , (4.4) 
(1 + Az) z 

with the free parameter z. One can combine both objectives for a special set of 
parameters (z = —1/16) 

16 2 1 u ^ 

oii = y. ai2 = "9' a22 = l44- (45) 

The corresponding action was introduced by Garcia Perez et al. M under the name 



Square Symanzik action. In ref. pjl the asymptotic ((3 — > 00) anisotropy coefficients 



of this action have actually been calculated. For 577(3) one obtains 

c CT (0) = 0.09829281 , c T (0) = -0.03008298 . (4.6) 

It is therefore obvious to investigate the properties of this action further and in par- 
ticular to determine its anisotropy coefficients non-perturbatively. For comparison 



19 



Action 


R D 


Ri 


r = Rn/Rj 


(2 x 2) 
Square Sym. 


1.088 
0.957 


0.99 
0.91 


1.099 
1.052 



Table 3: Derivative (Rd) and integral method (Ri) ratios at N T = 4 for the (2 x 2) 
and the Square Symanzik actions. 

we have as well considered the somewhat simpler (2 x 2)— action, which is improved, 
but its perturbative anisotropy coefficients are unknown. 

In order to apply the integral method we need again the high temperature limits 
of energy density and pressure. For the Square Symanzik action we have derived 
the expansions corresponding to eqs. ( |2.9|) and ( p.lOj ) 

« 1 - 0.1732 - 0.0457 ^) + O (iV^ 10 

Ri{Nt) = 1 _lfM'_»(M s '«(M ! + oM( , 

n ; 105 \N T J 693 \N T J 4729725 \N T J \ T J K 

I Q.07G [ - 0.115 ( - 0.089 ( ZX + O (N~ 10 



NJ \NJ \N T 

By chance, the ratio Rd(N t ) is even 0(a 4 ) improved. Since we intend to make use 
of Monte Carlo simulations on iV T = 4 lattices we list in Table |3| the respective 
numerical results for the complete 1/N T — expansions of the ratios. The numbers for 



the (2 x 2)— action have been taken from ref. [12]. 



4.1 Results for the (2 x 2)— action 



The thermodynamics of the SU(3) (2 x 2)— action has already been investigated 
in detail by simulations on 24 3 x 4 and 24 4 lattices in ref. |T2|] . We can therefore 
immediately apply the integral method to compute the anisotropy coefficients using 
the corresponding plaquette expectation values and the j3— function. For N T =A the 
critical coupling was found to be (3 C = 4.3995(2) . The behaviour of c a is similar 
to that shown in Fig. [I| for the Wilson action. With increasing (3 c a is slowly 
decreasing. In the neighbourhood of the critical point the integral method leads 
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p 


4.4 5.0 5.9 


C T 


0.226(08) 0.157(16) 0.115(16) 
-0.190(08) -0.096(16) -0.067(16) 


c a {k = 1) 
c r {k = 1) 


0.236(03) 0.195(07) 0.158(13) 
-0.200(03) -0.134(07) -0.110(13) 



Table 4: The matching results for c a and c T for the (2 x 2)— action. 



again to numerical difficulties. In Fig. [T3| we compare the integral method results 
with matching results, which were obtained on anisotropic 24 4 lattices with 2000- 
4000 Wilson loop measurements for (3 = 4.4 and 5.0 and 600 for (3 = 5.9 . We note, 
that due to the inclusion of (2 x 2)— Wilson loops in the action, less links can be 
integrated to increase the accuracy of the measurements. In Table ^ we present 
the values for c a and c r from the matching procedure. Again we observe, that the 
matching with fixed k = 1 leads to somewhat higher values. We find agreement 
between the different methods, but like for the case of the Wilson action only after 
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Figure 13: The anisotropy coefficient c a of the (2 x 2)— action from the integral (blue 
squares) and matching methods (circles). Results from fixed k = 1 matching fits are 
shown as green squares. 
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taking into account the cut-off correction factor r from Table ^. 



4.2 Results for the Square Symanzik action 

In order to determine the /3— function for this action we have performed simulations 
on 24 4 lattices and measured the plaquettes and the string tension as in ref. [2~0|. 
The lattice result fr(/3) and the physical string tension a are related by <r(/3) = aa 2 . 
The f3— function is then derived from 

_ / / — — ^ -i 

dg~ 2 a d/3 



da 2N r da 




(4.9) 



In Fig. [14] we show our measured values for Ja((3) together with a fit from which 
we have deduced the numerical values of the (3— function. On a 24 3 x 4 lattice we 
have then measured the generalized plaquettes P a and P T as well as the Polyakov 
loop and its susceptibility in the range 3.9 < (3 < 6.0 . On the average we took 
4000-8000, close to the critical point 18000 measurements. From the peak of the 
susceptibility we estimate the following critical coupling for N T = 4 

P c = 3.9820 ( If ) . (4.10) 

As a side product of our calculations we obtain the ratio 

-4= = \ = 0.643(4) , (4.11) 

° iV T ^(&) 



which is in agreement with results found for other improved actions Using 
the above mentioned data and r from Table ^ we have again applied the integral 
method to find the anisotropy coefficient c a . It is shown in Fig. |T3| together with 



the asymptotic value c CTi00 , eq. Q4.6D , and the results of our matching analysis. The 
latter are also listed in Table [5| They were obtained from Wilson loop measurements 
on 24 4 lattices. Like for the (2 x 2)— action less link integrations than in the Wilson 
action case can be done. Comparing the different results in Fig. [T5| we find again 



the same general (3— dependence as for the other actions. Also, a sizeable difference 
between the different matching options for k is found: the k = 1 results show an 
even stronger influence of the Wilson loop before, the k ^ 1 results agree 

nicely with the integral method data. We note that c a in the Square Symanzik case 
is about half as big as in the Wilson case. 
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Figure 14: The square root of the string tension, ycr(/3) for the SU(3) Square 
Symanzik action. The line is a fit to the data. 
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Figure 15: The anisotropy coefficient c a of the Square Symanzik action from the 
integral (blue squares) and matching methods (circles). Results from fixed k = 1 
matching fits are shown as green squares. The dashed black line shows c CTj00 . 
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p 


4.0 


4.5 


4.75 


5.0 


5.5 


N 

1 > meas 


2000 


800-1000 


700 


500-1500 


200 




0.149(11) 


0.143(17) 


0.132(20) 


0.119(18) 


0.109(21) 


C T 


-0.102(11) 


-0.065(17) 


-0.042(20) 


-0.031(18) 


-0.033(21) 


C a (k = 1) 


0.175(03) 


0.158(06) 


0.155(09) 


0.157(09) 


0.120(13) 


c T (k = 1) 


-0.128(03) 


-0.080(06) 


-0.065(09) 


-0.069(09) 


-0.044(13) 



Table 5: The matching results for c a and c T for the Square Symanzik action. N, 
is the number of Wilson loop measurements . 

5 Summary and conclusions 



One of the major objectives of our paper was the verification of the equivalence of the 
integral and derivative methods for the calculation of thermodynamic quantities like 
energy density and pressure. The major obstacle in the comparison of the methods is 
the cut-off dependence on finite lattices. In the high temperature limit the deviations 
of energy density and pressure from their respective continuum limits due to the cut- 
off effects can be calculated. They manifest themselves as a dependence on 1/iVj? in 
the thermodynamic limit, and they are different for the two methods. 

The anisotropy coefficients as such must be finite size independent. Their cal- 
culation allows then for a check on the internal consistency of the integral method 
and of its equivalence to the derivative method. Moreover, its non-perturbative 
behaviour is of importance for the approach of thermodynamic quantities to the 
continuum limit. 

By taking into account the major cut-off effects through the correction factor r 
we have shown for the Wilson action that one obtains similar results for N T =4, 6 
and 8 which agree, inside error bars, also with the matching results and with the 
results of Ejiri et al.||15||. The c CT — data of Klassen|TS[ are however definitely higher. 



This may be connected to the difference in our matching methods, the accuracy 
of the data and the lattice sizes used. We found it essential for the application of 
the matching procedure, that link integration of the Wilson loops was carried out, 
wherever this was possible, and that ratios instead of Wilson loops directly were used 
for the matching. Thus we avoided perimeter and corner dependencies. Remaining 
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effects of this kind are probably responsible for the difference of k = 1 to k 1 
matching fits. The fact that £(7) was linear near 7 = 1 helped in the determination 
of the derivative (d^/dj)^. 

We have repeated our analysis for two improved actions, the (2 x 2) and the 
Square Symanzik action, and find similar behaviours of the corresponding anisotropy 
coefficients. For all investigated actions the non-perturbative c a decrease with in- 
creasing (3 and approach their respective weak coupling limits from above; the results 
from the matching and integral methods are in all cases compatible with each other. 

In addition we determined for the Square Symanzik action the deconfinement 
transition point (3 C for iV r = 4 , the ratio T c j yfo and the (3— function in terms 
of the string tension. As an analytic result we derived the thermodynamic high 
temperature limits for this action. 
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